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ABSTRACT 

Tho  equations  of  notion  of  a body  with  slight  asymB  tries,  such  as 
a bonb  with  a bent  tail,  are  derived,  uxLng  linear  aerodynanlos. 

The  "resonant"  effect  of  coincidence  between  yaw  and  spin  rates  is 
discussed.  Conditions  that  the  spin  and  yaw  be  "looked  in"  (spin  rate 
equals  yaw  rate)  are  found,  applicable  in  the  stapler  cases. 


I.  INTRODUCTION 

Because  of  manufacturing  tolerances*  damage  in  handling,  etc.,  some 
projectiles  differ  appreciably  in  form  or  weight  distribution  from  the 
perfect  symmetry  assumed  in  ordinary  ballistic  theory.  This  report  sets 
up  the  differential  equations  of  motion  of  a projectile  with  certain  snail 
asymmetries,  and  discusses  qualitatively  some  features  of  the  notion.* 

To  fix  the  ideas,  the  analysis  will  be  made  for  a rigid  bomb  with  a 
symmetrical  body  whose  principal  axis  of  inertia  coincides  with  the  axis 
of  symmetry.  The  tail  will  be  assumed  weightless,  with  its  center  of 
pressure  not  on  the  axis  of  the  body  and  its  lines  of  sero  lift  and  moment 
not  parallel  to  this  axis.  If  the  body  is  unbalanced,  it  can  be  considered 
aerodynaraically  as  a balanced  body  with  an  appropriate  fictitious  tail, 
and  the  analysis  will  have  the  same  form. 

II.  SYSTEM  CF  FORCES  AND  MOMENTS 

Let  x^,  Xg,  x3  be  a right  handed  system  of  unit  vectors  on  axes 

rigidly  attached  to  the  projectile,  with  origin  at  the  center  of  gravity, 
x-  pointing  in  the  forward  direction  along  the  axis  of  symmetry  (hereafter 

called  the  axis),  and  the  other  axes  so  chosen  that  the  center  of  pressure 
of  the  tail  is  at  - tdx^  + hdxg,  where  d is  the  di  ameter  of  the  body.  Let  s 

be  a unit  vector  along  the  trajectory.  Let  ft  be  the  vector  angular  velocity 
of  the  projectile,  measured  with  respect  to  fixed  axes.  When  resolved  on 
the  moving  axes,  ft  £3*  bet  H be  the  vector 

velocity  of  the  center  of  gravity,  u - u *.  See  Figure  1. 

yW  hr 

,The  line  of  aero  normal  force  for  the  tail  will  be  the  unit  vector 
/l  Xl  ♦ /2  Xg  ♦ £3  where  i3  nearly  unity. 


* 

Somewhat  less  general  treatments  of  the  problem  have  been  given  before; 
see  Ref.  3. 


Center  of  pressure  for  normal  force.  Since  the  analysis  is  intended 
chiefly  to  suggest  the  form  of  the  system  of  forces  and  moments,  not 
numerical  relations  among  the  members,  the  fact  that  the  center  of 
pressure  for  Magnus  force,  say,  is  not  at  this  point  can  be  taken  care 
of  by  the  siae  of  the  corresponding  aerodynamic  force  coefficient.  It 
is  assumed  that  all  centers  of  pressure  for  the  tail  lie  on  a line 
through  the  center  of  gravity. 


"""By  the  normal  force  on  the  tail  will  be  meant  the  difference  between  the 
normal  force  on  the  complete  bomb  with  straight  tail  and  on  the  body  alone, 
as  measured,  for  example,  in  a wind  tunnelo 
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The  following  forces  and  moments  will  be  assumed  to  act: 

1.  Gravity. 

2.  (*)  A foroe  on  the  body  whose  axial  and  cross  components  are 
respectively  given  by  th9  axial  and  cross  components  of  the  velocity  of 

a representative  point  on  the  axis;  i.e,,  the  components  of  uz  ♦ L Q x x^, 

where  L is  a constant  length.  This  force  includes  axial  drag,  normal 
force,  and  cross-spin  (pitching)  force. 

(b)  A corresponding  force  on  the  tail,  with  the  resolution 
along  and  perpendicular  to  the  line  of  zero  normal  force.  The  velocity 
components  are  found  by  writing  the  velocity  of  the  center  of  pressure 
of  the  tail  ass 

u(s  • »’)  jj'  ♦ u £z  - (z  • »•)  s*J  ♦ Q x (-tdxj  ♦ hdXg). 

3.  (a)  A force  on  the  body  given  by  0 x (u  i * L O x x^),  the 

Magnus  force  and  Magnus  cross  force.  This  generalizes  the  Magnus  force 
on  a baseball,  Q x u. 

(b)  A similar  force  on  the  tail,  given  by  Q x u(z  • z*)  z * 

♦ u £(*•(*•  *')  *'  J ♦ Q x {-tax^  ♦ hcbcg)^  • 

U*  The  moments  of  these  forces,  given  by  x^  x (u  z + L Q x x^),  etc. 

5.  Other  moment ss  Spin  damping,  spin  accelerating  if  the  fins  are 
canted,  and  moments  resulting  from  non-coincidence  of  the  lines  of  zero 
normal  force  and  moment. of  the  tail. 

To  compute  uhe  detailed  form  of  these  forces  and  moments,  write: 

uz-u^x^  + UgXg  ♦ u^  Then  Uj/u  ■ cos  6,  where  6 is  the  angle  of 

2 2 2 2 

yaw,  and  (ug  ♦ u^  )/u  - sin  6.  Then  2(a)  gives  farces  proportional  so 

^ and  (ug  =;>  ♦ u^  Xj)  ♦ L ( D -j  Xg  - & 2 ♦ 

—1 

The  factors  of  proportionality  must  have  the  dimensions  MT  , and  aero- 
dynamic theory  and  experience  indicate  that  they  should  be  pd2  uK,  where 


p is  the  density  of  the  sir  and  K is  a dimensionless  constant.  Then 
these  forces  will  be  written 

(1)  - pd2  mu  Kjjja  x^  - pd2  uK^Cug  Xg  ♦ u3  x3)  - pd3  uK^  (U3  Xg  -ttg  x3), 

the  axial  drag,  normal  foroe,  and  crosa-apin  force,  with  the  subscript 
one  indicating  that  ^he  body  is  referred  to.  The  minus  signs  make  IU., 
and  positive. 

In  working  out  2(b)  it  will  be  assumed  that  Ug/u,  u^/u,  h>  and  f3 

are  so  smell  that  the  product  of  two  of  them  is  negligible.  Also,  /, 

2 2 ± 

will  be  taken  as  unity.  The  coefficients  - pd  u K^,  -pd  u and 
pd2  u will  be  applied  to  u(a  . a1)  a*,  u [ a - (a  . a')  a1  “J  , and 
Q x (-tdx^  + hdxg)  respectively.  The  result  is 

-pd2  UU!  (*1  ♦ t2  58  * h ^ ^ *N2  (U2  £2  + U3  & 

(2) 

♦ pd2  ui^  ^ (/2  Xg  ♦ x3)  - pd3u  Kgg  ( n 3 3^  - ft2  x3  ) 

* pd3  u Kg2  (h/t)  ( ^ 3 ^ i £3^* 

The  Magnus  forces  in  3(a)  are  proportional  toi 

Q1  ^v3  £2  " "2  £3)*  ( ^'3  £3)*  ( a2  ^2  + ^3-3)# 

2 2 

and  ( ft  ? ♦ ) Xy  Tine  first  gives  the  ordinary  Magnus  force, 

pd  Q.  ^ (^  Xg  - Ug  x3).  The  second  has  the  same  form  as  the  cross-spin 

foroe,  except  for  the  factor  instead  of  u.  Since  the  assumption 

cos  6-1  will  be  made  later,  this  force  will  be  included  in  the  cross- 
spin  force  of  2(a).  The  third  term  gives  the  Magnus  force  due  to  crostr 

apin,  pd1*  jjn  O x ( Q g Xg  ♦ 

li  2 2 

order  in  small  quantities}  it  would  give  a farce  pd4  K ( 2 ♦ Cl  3 ) 

j.  2 2 2 2 2 

or-  pa*  u K ( Cl  g /u  + ^ ^ /“  )*  where  the  parenthesis  is  of  the  order 
*2 

of  6 , so  it  will  be  dropped.  Then  the  Magnus  forces  on  the  body  will 
be  written 

(3)  - pd  O'  ^ (u3  Xg  - Ug  x3)  ♦ pd  ^ ^ 2 £2  + ^ 3 —3^ * 


H 3 x3).  The  fourth  term  is  of  second 


7 


In  computing  the  forces  from  3(b),  the  velocities  along  and  normal 
to  the  vector  £1  * / 2 £2  + 5,3  wi^1  k«pt  separate.  Thus  (^  . z1)  z' 

will  be  written  x^/u,  neglecting  small  terms,  and  z - (z  «,  £')  js!  will 

be  written  (ug/u  - /g/u)  Xg  ♦ (u^/u  - /y\i)  Xy  the  term  - (Ug  j^/u 

♦ jfy'u)  Xj_  being  dropped  because  it  is  of  second  order  in  small  quantities. 

Then  3(b)  gives  terms 

^1  ( n 3 % ~ ft  2 £3^ * “ ^ 1 (^3  “ u2  £3)*  ^ 1 U1  (^3  52  “ ^2  £3^* 


-td  ^ ( O 2 ^ 3 £3)*  end  “ td  (h/t)  & Xgj  terms  in  0 g2, 

ft  ^2,  ar>d  0 g (h/t)  have  been  dropped  beoause  (h/t)  is  regarded  as  3m all. 

The  first  term  will  be  absorbed  in  the  K term  as  before,  the  second  gives 

the  ordinary  Magnus  force  on  the  tail,  the  third  gives  the  Magnus  force  due 
to  the  tail's  line  of  zero  lift  not  being  parallel  to  the  axis,  the  fourth 
is  the  Magnus  cross-force,  and  the  fifth  is  a Magnus  force  resulting  from 
the  fact  that  the  tail  has  a velocity  hd  across  the  air  stream  when 

the  missile  is  spun.  Then  the  Magnus  forces  on  the  tail  will  be  written 

- pd3  ft  ^ Ky2  (uj  Xg  - Ug  x^)  - pd^  ft  Kjpg  ( ft  2 Xg  + ft.  3 £3) 


(M 


pd3  ux  Cl1  Xpg  (/3  Xg  - /g  x3)  - pdh  (h/t)  tKjyg  ft  x2  Xg. 


In  adding  (1),  (2).  (3),  (U)«  term!;  of  the  same  form  will  be  combined 
and  a coefficient  X with  no  numerical  subscript  will  be  assigned.  It  is 
not  intended  to  imply,  for  exampls,  that  K_,  below  is  Kp^  ♦ Kpj  but  merely 

that  « force  proportional  to  ft  ^ (u^  Xg  - Ug  x^)  may  bo  expected  to  appear. 

(Kpi  ♦ Kpg  may  be  a useful  first  estimate.)  Then  the  aerodynamic  for C6  on 

the  projectile  will  be  assumed  to  be* 


2 2 3 

- pd  uu1  XDA  - pd  uKjj  (ug  Xg  + u^  x^)  + pdJ  uKg  (ft^  -flg  x^ 

- pd3  Kp  ft  1 (u^  Xg  - Ug  x^)  + pd11  Kjp  0 1 ( ft  g Xg  ♦ x-j) 
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(5) 


+ pd  uu^  K|2  (/jj  2g  * «3  23^  + ui  ^1  ^2  ^3  22  ” ^2  £3) 

~ pd^  P^Kjyj.2  ^*1  2g  + pd  u^ti2  ^ ^ 3 2i  * ® i £3)* 


Hero  p - h/t  and  - lt^  - K^,  the  lift  coefficient  of  the  tail. 

The  moments  of  the  farces  about  the  center  of  gravity  can  be  fcund  thus: 
The  farces  (1)  and  (3) , anting  on  the  bidy,  are  assumed  to  act  through  points 
on  the  axis.  Then  their  moments  can  ts  found  by  taking  the  vector  products 
with  ju,  raising  the  powers  of  d by  one,  and  changing  the  aerodynamic  coef- 
ficients to  allow  for  the  positions  of  the  centers  of  pressure.  If  the  new 
K's  are  given  the  additional  subscript  M,  the  result  is: 

<j  J 

pd  ^MNl  ^u3  22  * 23)  “ Pd  “^MSl  ( ^ j 2g  + ^3  £3) 

(6) 


Pd1*^! 


Q 


1 (Ug  ^ ♦ u3  x3) 


pd?  K 


MFl 


ftl< 


ft 


3 % 


- ft 2 x3). 


The  forces  (2)  and  (U)  are  assumed  to  act  through  points  on  the  vector 
- td  x-^  ♦ hd  Xg.  Therefore,  this  vector  must  be  used  instead  of  to  find 

the  moments.  The  result  is: 

- pd3  utKjg  (v^  Xg  - u2  Xj)  ♦ pd^  utKgg  ( Q 2 2g  + ^3  23) 

- pd  G ^ tKj<2  (u2  £g  + u3  £3)  + pd  0 "*XF2  ^ ^ 3 22  **  ^2  2«) 

♦ pd"  ur^  tKTp  (/j  Xg  - Xj)  - pd  ft^  tKp2  ^2  22  + ^3 


♦ frxl3  | (-utKffi  u3  ♦ uuj^  tKj^  /3  - dutK^  ftg  ♦ d ft  3 tKp2  u2 

“ d ft  ^ tXjrp2  ^ 3 “ du^  0 ^ tKp2  /j*  ” Pdu  0 ^ tKg2)  x^  “ du  ft  ^ tK^  x 2 

* (UU1  tKDA2  +d2t^  l2  W 23}  * 


| 

i 
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Some  small  terms,  e.g«,  in  £ Kgg  and  pK^g  /g,  have  been  kept;  not 
because  they  are  effective  but  because  they  are  easily  seen  plysicolly. 

When  these  moments  are  added  to  those  in  (6),  corresponding  terms 
will  be  lumped  together;  e.g.,  pd^  u(Kj*,^  - tKjjg)  Xg  - Ug  x^)  will 

be  written  pd3  uKM  (u^  - Ug  x3).  The  notation  will  be  according  to 

the  following  scheme; 

Force  Coeff.  Associated  Moment  Coeff. 


K 


N 


K 


M 


* 

hr 


h 

h 

*XT 


The  coefficients  of  moments  due  to  asyanatry,  (those  containing  fat  fat 
or  h/t)  will  bs  left  as  they  are;  not  because  tKgg,  for  example,  is 

neoessarily  the  correct  coefficient,  but  to  show  the  derivation  and 
help  in  making  rough  estimates,  The  system  of  moments,  given  ty  U. 
above,  is  them 

pd3  | uXM  (u3  Xg  - «o  x,)  - duKH  ( O 2 Xg  ♦ n3  x3) 

“■  d O ^ K^,  (tig  Xg  + u3  f 3)  — d Q ^ Kjj  ( Q 3 Xg  * 0 g 

♦ uu^  tKj^g  ( fa  £2  " —y'  “ ^U1  l ^ ig  * fa  —3^ 

(7) 

♦ Ppd  (uux  tKT  9 “ utKjjg  u3  * dutKgg  0 g * d ft  ^ tKjpg  Ug 

- d H 1 tK^pg  O 3 - duj.  ii,  tKpg  /g)  3^ 


- du  n1  tKgg  xg  ♦ (u^  tKjj^g  ♦ d2t2r*  x2  k^)  x3  f . 

In  addition,  there  will  be  assumed;  (a)  Spin  damping  moments, 
proportional  to  H ^ x^  for  the  body  and  £ 0 . (-td  x^  + hd  Xg)"^ 

(-td  v-^  + hd  Xg)  for  the  tail,  the  aerodynamic  coefficients  being  called 
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KA1  KA2* 
[ ^ - (r/u) 


(b)  A spin  generating  moment  on  the  tail*  proportional  to: 
[n.  (-tdx^  + hdxg)]  | (-tdx^  + hdXg),  where  fc  is  the  angle 


of  cant  of  the  fins  and  r is  the  distance  from  the  axis  of  the  tail  to  the 
center  of  pressure  of  one  finj  the  coefficient  will  be  Kg,  (c)  If  the 

axis  of  zero  «c»ent  of  the  tsil  does  not  coincide  with  the  axis  of  zero 

3 2 

lift,  an  additional  moment  pd  u Kc  *1  * °2  * n3  *3^  wil^  be 

2 2 

assumed,  where  rig  ♦ n^  "1  approximately  and  n^  is  small.  When  second 
order  terms  are  dropped,  these  moments  are : 


pd3  u I (uKg  € - **%  A x - dKA  a ! ♦ uKc  1^) 

(8) 

♦ (-upKg  < + rpKjj  ♦ dpKA  Ci1  ♦ uK,,  r^)  Xg  ♦ uKc  n3  x. } ] . 


The  second  and  third  terms  in  each  parenthesis  will  combine,  as  will 
the  first  and  fourth,  so  only  the  first  and  second  will  be  written.  In 
say  actual  problem  these  combinations  must  not  be  forgotten.  Then  (7) 
and  (6)  give  the  moments  considered  hereafter. 

The  force  of  gravity  will  be  represented  by  the  vector 
(9) 


£ “ 81  5i  + % *2  4 83  £3’ 

III.  EQUATIONS  OF  MOTION 

These  will  bet 

m (d/dt)  (u^  x^  + Ug  Xg  ♦ x^)  ■ stun  of  the  forces,  and 

(d/dt)  (Afl^x^  + BOgXg  ♦ B Cl  ^ x^)  “ sum  of  moments  about 
center  of  gravity, 

where  m is  the  mass,  A is  the  axial  moment  of  inertia,  and  B the  moment 
of  inertia  about  a transverse  axis  through  the  center  of  gravity.  Befo:;e 
writing  down  the  equations  in  component  form,  the  force  equations  will  be 
divided  by  m and  the  moment  equations  by  B.  Wherever  B divides  an  aero- 
dynamic coefficient  it  will  be  replaced  by  mk^d^  where  k is  the  transverse 
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radius  of  gyration  in  cslibarsj  and  ary  combination  (pdV«)K  will  bo 
raplaced  by  Ji  Tha  equations  of  motion  ara  than,  using  (5)»  (7),  (8) 
and  (?), 

(a)  ^ - tx2  A3  ♦ Uj  0 2 - - (J^/d)  us^  ♦ upJ^  0 ^ ♦ gx, 

(b)  + - “u(Jjj/d)  Ug  - uJ£  ^ - Jj,  O 3 u^ 

♦ dJjF  n i a,  * ,JUl  (jL?/d)  /2  * «!  o 3.  /3  - pvujj.jt)  x*  • gj, 

(c)  ^3  “ \ & 2 * ^ ^ 1 " “ u ( Ji/d)  u3  + uJs  ^2  + ^ 1 ^ 

+ °1  °3  * UU1  (JL2/d^  h “ *1  ^2  fll  *2  " uPJS2ftl+ 


(d)  H x - (B/Ak2d2)  ^ u2  J£  i - udJA  D1  + p [ tJ^  /3  - utJN2  u^ 

- utdJ^/)  2 * tdJF2  “2  ” ^ 1 JXF.2  ^ 3 “ U1  ^ 1 ^*^2  ^2 

- Pu  0 ^ 1 ^ » 

(e)  O 2 - ^i  ^3  (1  - A/k)  ■ '^[uJM  u3  " 2 " ^ 1 dJT  “2 

^ 1 d + uu^  *3  " ul  ^ l tdJpg  ^2  ” u 0^e6 

♦ udpjAa1-pua1  tdJS2]  /k2d2, 

(f)  03  + fjg  (1  - A/B)  ■ f-  uJjj  Ug  - udJH  A ^ dJ,j,  U2 

2 2 

^ ^ d Jjt  2 **  ^ fl  ^ “•Vj  h * “ Kc  ”3 


♦ B (u^  tJDA2  ♦ Clf  d2  t2  Jjp2)]  A2  d2 


To  these  must  be  added  the  equations  for  g,,  g~,  g-,  obtained  by 
differentiating  ^ 


gl5i  + g222  + g3r3“  constant* 

(g)  - gj  A 3 * g3  ^ 2 " °» 

(10) 

(h)  gg  - g~  i)  ^ ♦ g^  O 3 ■ 0, 

(i)  83  ” 8X  a 2 + *2  0 1 " 0i 
and  the  equation 

(10) 

(j)  u,2  ♦ u22  ♦ u,2  • u2. 

Equations  (10)  give  the  motion  of  the  projectile  under  the  assumptions 

made. 


The  algebraic  signs  have  been  chosen  thus: 

KK’  KH.j  l?»  KXF  *re  positive  if  the  corresponding  forces  act  ahead  of  the 

center  of  gravity}  e.g.,  is  negative  for  a statically  stable  bomb. 
Aerodynamic  coefficients  for  the  tail  alone  -are  positive,  provided  the 
Magnus  force  is  given  in  direction  by  (3pin  vector)  x (velocity  vector), 
and  acts  behind  the  under  of  gravity. 

Before  proceeding  with  the  analysis,  dimensionless  variables  will 
be  introduced.  Let: 


v - u/d,  vi  - u^d,  i - 1,  2,  3, 


CO.  = A)  ±/\, 


i «*  1,  2,  3,  and  a 


o 


vdt. 


Since  the  system  of  aerodynamic  forces  and  momsnts  ia  only  plausible 
for  small  yaws,  let  uu/u  " cos  6 ■ 1.  Let  primes  denote  differentiation 

vitn  respect  to  a.  Then  the  equations  of  motion  become  i 

(a)  y^/y  - (^2  03  - v3  oo^/v  - - JDA  ♦ pJs2  oxj  ♦ Qx, 

(b)  v2'/7  " (y3  “1  " Yl  “3VV  - - JN  Vj/v  - Js  - Jp  “1  V3A 

2 

+ *1  ®2  * ^L2  * ^F2  ^ ^3  ” mi  + ®2* 

(11) 

(c)  v3'/v  - <*2  “ v2  <\)A  ■ “ JN  v3/v  JS  “1  v2^v 

+ JXF  ^ “3  + JL2  h “ "^2  ®1  ^2  “ ^S2  “1  * °3* 

(d)  o^'  +6 ^ v'/v  ■ (B/Ak2)  £ Jg  « • + ptU^  £3  “ JN2  v3/v 

- JS2  «2  ♦ «i  JF2  v2/v  • ^1  Jff2  cij  - l2  - JJc^  JS2)  ] , 

(e)  c o^'  + v‘/v  - 0^  toj  (1  - A/B)  « [ JM  v^/v  - JH  - 03^  JT  Vg/v 

“ *1  JXT  ®3  * tJL2  h “ ®1  tJF2  ^2’^E£  +PJAa>l"pCftl  tJS2  ] ^k2* 

(f)  O^'  ♦ CO3  Y'/v  ♦ cou  (1-A/fe)  - [ - JM  Y 2/v  - JR  «3  - JT  Y3/y 

♦ ®1  Jjt  tJF2  /3  ♦ Jc  «3  ♦ ptJDA2  * P^2  t2  j^]  A2; 

(g)  Q^»  ♦ 2Qj^  v'/v  “ Q2  “3  + ^ ^ • 0, 

(h)  Qg*  -v  2Q2  y*/y  - Q3  0^  ♦ Ql  - 0, 
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(i)  G^'  * ?G^  v'/v  - Gj.  ci^  + Gg  « 0; 

(J)  v2  - v-j2  ♦ v22  ♦ v32, 

(k)  w'  - v^  ♦ v2  v'2  + v'3 

Now  velocities  are  in  calibers/sec.,  angular  velocities  in  radians/ 
caliber  of  travel,  and  derivatives  are  taken  with  respect  to  calibers  of 

2 2 2 

travel;  and  (v2/v)  » (v-,/v)  - sin  6. 

IV.  THE  YAWING  MOTION 

Let  (v2  ♦ ivy/v  * - £ , £iv>  * ica^  ■ p,  and  Gg  + i-G^  - T . Then  from 

(11)  (b)  and  (c),  (a)  and  (f),  and  (h)  and  (i), 

(a)  £ ' + i (v'/v-  ♦ ioij)  ♦ in(v^/v)  - £ (-JN  ♦ ia^Jp)  - |i(Us  + o^J^) 

" + ^3^  + Pt,JxF2a>l  “ ^ ? 

(12) 

(b)  p'  ♦ pv'/v  ♦ in^U-A/B)  - 4 (iJK  ♦ - ^(JH“icaI  Jxt^ 

- (/2  ♦ 1 ♦ c^tJyg)  ♦ U0n3  ♦ .s(itJDA2  ♦ i^2^JXF2 

“ JE€  * JA®1  " “ltJS2^» 

(c)  r • ♦ 2 r v'/v  ♦ icojT  - iO^p  ■ 0 

Then  the  complex  numbers  ^ , n,T  represent  the  sine  of  the  angle  of  yaw, 
the  component  of  the  angular  velocity  perpendicular  to  the  x^  axis  (cross 
component),  and  the  component  of  gravity  perpendicular  to  the  x-^  axis, 
respectively. 


1$ 
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: 
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l 
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It  is  now  convenient  to  change  to  axes  not  fixed  in  the  projectile* 
The  substitutions 


(13) 


t • £ exp  (i 
G - P exp  (i 


l k exp  (i 


rotate  the  axes  back  as  fast  as  the  projectile  spins  forward.  Then  by 
proper  choice  of  the  aero  point  of  a,  the  axes  can  be  taken  as  (say) 
forward  along  the  projectile's  aria,  and  approximately  down  in  a vertical 
plans  containing  the  axis,  and  to  the  left,  horizontally,  far  small  yaws 
and  over  moderately  long  pieces  of  the  trajectory* 

Equations  (12)  become 


(a)  f ' ♦ / v«/v  ♦ ipvj/v  - £(-Jn  ♦ io^Jp)  - p(iJg  ♦ “ G 

a 

■*  £ ( J| £ “ i°^L  "*  i^l^S2  " 1 exP  (1  P ^da), 

Jo 

(1U) 

(b)  p'  ♦ p(v'/v  - ia^A/B)  - | c (iJM  ♦ c^Jj)  - p(JH  - icc^  J^) 


♦ [ - (/2  ♦ i/jKitJj^  ♦ eo^t^)  *iJcn3 

a 

♦ «i<-JDA2  * * ¥*  JA“l  - “ltJS2>  ] eIP(i  J 

(c)  G«  ♦ SGv'/v-iC^p  - 0. 


To  get  the  differential  equation  for  the  yawing  motion,  differentiate 
(lU)(a),  substitute  for  p'  and  O'  from  (ib)(b)  and  (c),  and  then  substitute 
far  p from  (lU)(a).  Products  of  J's  will  be  omitted  from  coefficients  of 
g and  g ' in  the  final  equation,  but  retained  (to  the  second  order)  in  the 
right  member.  The  exponential  terms  in  (lb) (a)  and  (b)  will  be  denoted  by 
and  qg  for  convenience.  Since  v^/v  ■ u^/u  ■ cos  6,  it  will  be  set  equal 


I 

i 

| 


! 

i 
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to  unity.  The  coefficient  of  p in  (lU) (a)  Kill  be  taken  as  constant,  since 
(Jj.  *j)  1 can  hardly  be  appreciable  in  most  cases.  The  expressions  q^  and 

will  be  treated  as  of  second  order.  It  is  shown  in  Ref.  2 that  good 
approximations  to  <3^  and  0 are 

* O 

0^  * - (g/vM)  sin 


0 » (g/v2d)  sin  0 £ * (g/v2d)  ccs  S, 


where  © is  the  angle  from  the  horisontal  to  the  trajectory,  positive  up-_ 
wards;  these  approximations  will  be  used  in (14)  (c)  to  compute  Q*  and  (g/trd) 
will  be  called  y and  as  stated  small.  All  J's  will  be  treated  as  constants; 
their  derivative"  can  hardly  be  larger  than  second  order,  except  in  very 
unusual  cases.  The  quotient  v'/v  will  be  snail  far  shell  and  bombs,  but 
may  be  fairly  large  for  rockets,  so  it  will  be  left  in  wherever  it  appears. 
The  equation  for  the  yawing  notion  then  becomes 

{*+£'[  2r'A  + J„  ♦ Jj/k2  ♦ Y sin©  - ic^A/B  - ico^Jp*  Jjy/k2)] 


♦ £ [-  Vk2  ♦ O'A)*  ♦ (v'/v)2  ♦ (v‘/v)(JN  ♦ JH/k2)  - o^2  (A/BK 


(15) 

- ic6^(A^)(v*/v  ♦ JN  ♦ Y ®in®)  ♦ io&^J^/k2  - ixo^ (▼•/*)  Jy] 


, O , p 

■ y cos©  (v'/v  • Jg/k  “ Y ®in®  ♦ io&^A/B  ♦ iJj^/k  ) 


♦ ^(v'/v  - ia^A/B)  ♦ - iqg  . 

(The  reader  may  note  that  in  the  simple  case  of  vertical  fire  in  a 
vacuum,  this  equation  is  incorrect,  because  the  term  y2  sin?©  in  the 
coefficient  of  £ has  been  dropped.  If  this  term  is  put  in,  .the  equation 
will  give  the  correct  result.) 


This  is  not  really  necessary  here;  the  important  assumption  is  that  (v,/v)  • 
is  negligible.  However,  in  combining  certain  force  terms  earlier  this  J 
assumption  was  made,  and  certainly  for  the  motions  of  greatest  interest  the 
yaw  is  small. 
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Except  for  some  changes  in  notation,  and  the  addition  of  the  terms 
in  q.,  and  qg,  this  equation  is  identical  with  equation  (1*23)  of  Ref*  1, 

as  corrected  in  Ref*  2* 

The  yawing  motion  will  have  (usually)  a transient  oscillatory  part, 
compounded  of  two  damped  circular  motions  with  rates 


(16)  (Affl1/2B)(1  ♦ c), 

where  ■ 1 - l/s,  s ■ and  l/s  is  < 1;  and  a "steady" 

yaw,  a particular  integral  of  equation  (15) » The  formulas  (16),  and 
conditions  for  the  damping  of  the  transient  motion,  are  discussed  in  Ref.  1. 

V.  RESONANCE 

The  steady  yaw  ran  be  broken  into  two  parts.  The  first  is  computed 
by  putting  q.  ■ q^  • 0,  and  is  shown  in  Ref.  1 to  give  a small  yaw  causing 

drift.  The  second  can  be  computed  by  putting  y a 0 in  the  right  member 
of  (15).  In  general,  no  convenient  analytical  expression  is  available  far 
this  result.  However,  it  is  useful  to  consider  the  case  of  a projectile 
moving  at  nearly  constant  velocity,  with  ca^  changing  slowly. 

If  y Id  set  equal  to  aero  in  the  right  member  of  (15),  the  equation 
can  be  written 


(17)  ♦ (C  ♦ iD)  £ • + (E  ♦ IF)  / - Q exp(  f ioojda). 

If  a solution  is  sought  in  the  form  exp(  f rda),  where  r is  a slowly 
changing  function,  substitution  gives.  J 


£m  Q exp  ( J'  ico^da)/  [ r'  ♦ r2  ♦ (C  ♦ iD)  r + (E  + iF)  j . 


Now  if  the  variation  in  Q,  C,  D,  E and  F can  be  neglected,  this  is  a 
solution  of  (17)  if  r - ito^ . Bui  the  denominator  is  the  left  member  of  • 
the  Riccati  equation  obtained  when  seeking  the  expression  for  the  transient 
oscillatory  yaw.  Thus  if  co.  is  one  of  the  natural  rates  (16),  the  denomi- 
nator is  very  small;  the  phenomenon  is  similar  to  the  resonance  of  a damped 


"if  ■>.  f oil  i;-;  fired  -t  largo  angles  of  elevation,  the  yaw  near  the  summit 
of  U:l:  trajectory  riay  bo  too  large  for  the  simple  theory  of  this  report. 


18 


oscillating  system  in  one  degree  of  freedom;  the  magnitude  of  £ is  re- 
stricted only  by  the  amount  of  damping.  This  resonance  can  be  illustrated 


by  plotting  | §/  [ - <a£  * ia^C  ♦ iD)  ♦ E * IF  ] | as  a function  of  eo^, 
assuming  that  is  negligible. 


Graphs  to  illustrate  resonance  in  the  motion  of  a fin-stabilised 
missile  appear  in  Figures  3 and  Ip.  To  avoid  a multiplicity  of  cases,  several 
simplifying  restrictions  have  been  made  in  preparing  these  grqpho.  The 
procedure  was  this:  v'/v  and  «\j  were  assumed  negligible,  and  qj  was 

written  ioo^q^.  Then  1 


Q - - [ t - ic^d-A/B)  ] (/2  ♦ i/3)  ♦ JQn3  ♦ fi  [tJDA2 

♦ i(JE6-  «a1JA)]e  [ t - i^(l-A/B)]  [iJF2(/2  ♦ i/3) 

* f<1JS2  * % * “1%.  say,. 

Now  a fin-stabilized  projectile  has  J < 0;  hence,  by  the  definitions 
after  (16),  o is  > 1*  Then  oft,  can  be  equal  to  only  one  natural  rate,  that 

* 2 2 

of  the  same  sign.*  If  this  occurs,  (16)  shows  m (-J^k  )/(l-A/b),  which 
for  ordinary  missiles  is  of  the  order  of  10”^.  The  values  of  of  most 
interest  will  then  be  near  .01,  so  the  term  in  ^2  waa  droPPc<^* 

Also,  if  is  very  small,  J£6  - oft^,  the  principal  part  of  ooj,  (cf.  (ll)(d)) 
will  be  very  small,  and  so  can  be  omitted  from  Q^.  Since  t is  likely  to  be 
2 or  more,  t - ito^(l-A/fe)  was  replaced  by  t.  Then  the  approximations 

«1  - - tJL2</2  ♦ l/j)  ♦ ■Jj'H  ♦ P«Dt> 

«5  * “(JF2  ♦ Pa!>. 


are  reasonable.  These  approximations  were  not  needed  in  making  the  graphs, 
hut  the  argument  above  shows  that  useful  results  can  be  obtained  near  the 
resonant  frequency  by  making  the  graphs  as  described  below. 


*If  the  projectile  is  spin-stabilized  with  > 0,  o>,  cannot  be  equal  to 
either  rate. 


1 9 


The  effect  of  is  shown  conveniently  by  plotting  the  ratio  of 
IV  [-  a^2  ♦ italic  * ID)  ♦ (E  ♦ IF)]  | 
to  its  value  when  <a^  - 0.  This  can  be  written  as 

|(B  ♦ it^CC  + iD)  eE  ♦ iF]|  . 

Now  if  the  rate  (Aaj^/2B)(l-*o)  is  set  equal  to  N,  it  follows  that 
- Jj/k2  - N2 

If  J^(l-A/B)  ♦ Jy/k  ♦ Jj/k  is  put  equal  to  X,  the  ratio  can  be  written 

I/2 

[l  - (ttjAXVb)]  /f  [l  - (aj/N)(A/fe)  - (o^A'/^Cl-A/B)]2  > (®LA)2(X/N)2}  , 

which  is  plotted  in  Figure  3 as  * function  of  »jA»  (XA  i8  replaced  by 
X/N  at  co^  » K.)  This  graph  indicates  how  the  effect  of  some  eccentricities 
is  amplified  by  spin. 

The  effect  of  Qg  can  be  shown  by  plotting 

(ttjAOAaon*  denominator  as  before), 

as  is  done  in  Figure  U.  The  size  of  the  yaw  produced  can  be  found  by 
multiplying  the  ordinate  by  j Qgj  j i.c.,  by  |t(Jpg  * | • 

VI.  LOCKING  IN 


A tendency  for  the  spin  and  yaw  to  lv  ck  ii.  ; i.e.,  for  the  projectile 
to  spin  at  the  resonant  rate,  so  that  the  angle  between  the  plane  of  yaw 
and  the  plane  through  the  axis  and  the  vector  hdxg  is  constant,  can  be 
shown  in  certain  rather  idealised  cases. 


Consider  a projectile  on  a straight  trajectory,  v*/v  constant,  with  a 
circular  yawing  nation  of  constant  rate.  (Far  small  spins,  Lhe  rates  (16) 
are  not  very  sensitive  to  changes  in  the  spin  rate.)  Figure  2 shows  the 
Eulerian  angles  6,  0,  ljrt  and  the  axes  X^,  Z-,  ly  fixed  in  direction  in 

opaoe,  with  origin  at  the  c.g,  and  X^  along  the  trajectory.  If  6 is  small, 

- $ cos  6 + ijf  approximately,  so  0^  ■ ♦ ifr* . From  Figure  2 it 

can  be  seen  that  the  complex  yaw  C ecu  be  written 


t ~ tci{^“  "/2) 


- X 
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and  further  that 

v2/v  ■ 6 sin  iff,  Vj/'v  - 6 cos  ijr . 

Now  it  can  be  seen  that  the  component  of  angular  velocity  perpendicular 
to  the  axis  is 

if*  ie^  ( 6 ♦ i6$0 . 

Since  b m 0 here,  the  component  wanted  is 
- ($(cos£  * sin0), 

from  which 

e4j  “ 6 0*  cos  y,  co^  ■ - 6 ff'  sin  y. 

Substitution  for  vg,  mid  the  ca's  in  (ll) (d),  and  neglect  of  0",  give 

(18)  y - y (a,  ♦ a2  sin  y)  + a^  sin  y ♦ a^  coat/r*  a^, 

where  ^ - - v'/v  - (*J'A  * hJF2jf2  + 

*2  “ (hb/k^CJ^  ♦ 0'J^). 

«3  “0'«2* 

au  - - (ho/ka2)(JM2  *JZl'Jg2), 

«5  “ - *•(▼'/▼)  ♦ [ JE<  - JA  ♦ h(JI2jT3  - %2l2  - PJS2)  ] A,2, 

2 ° 

where  ka  « B/Akc;  is  the  axial  radius  of  gyration  in  calibers. 
Equation  (18)  has  the  solution 


(19) 

a^  sin  f + 8} 

cos  tlr  + a-  ■ 0, 
1 5 

, r 2 2-172 

or 

sin  ( yr  * r)  ) 

- - a^/  [ a,2  4 au  J , 

where 

sin  Tj  ■ a / 

r2  2-j3^2  /r2  2 1 ^ 

l a3  + % J 'M5rvi‘)  * 

provided  that  the  inequality 

2s  2 2 

(20)  a^  « a3  + a^ 
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holds.  This  inequality  is  then  necessary  and  sufficient  that  a solution 
¥ ■ constant  exist}  i.e.,  that  a "locked-in"  state  be  possible.  Then  for 
h6  sufficiently  small,  locking  in  will  not  occur. 

It  is  convenient  to  define  the  angle  X by  the  equations  sin  X • 

2 2 ^ 2 2 2 ^ 2 2 ^ 

) , cos  X - (a3  ♦ a}i  - a^  ) /(a3  + a^  ) . Then 

(19)  is  satisfied  by 

f • - fj  ~ X end  f • - TJ  * X - n. 

l/2 

If  0*  is  small,  as  is  often  the  case,  (a^2  * a^2)  » h6 

approximately.  Then  sin  tf  ■ - 1,  *)  - - n/2  approximately.  Thus 
¥ • n/2  -X 

or  ¥ m * -W/2, 

or  from  Figure  2,  the  plane  of  the  bent  tail  differs  from  the  plane  of 
yaw  by  the  angle  X • 

Stability  of  the  Lockod-In  State.  Let  ¥c  he  a constant  solution 
of  (1C),  and  let  ~fr  • '{'Q  * T be  another  solution,  with  z small.  Then  if 
the  term  in  T r ' is  neglected,  (19)  gives 

t*  - ♦ *2  sin  fQ)  - r ( a^ cos  - a^sin  p0)  - 0. 

2 

The  discriminant  of  the  characteristic  equation  is  (a^  + aj  sin  ¥t ) 

U (a^cos  - a^sin  ^Q).  Then  for  stability;  i.e.,  for  r to  danp  out, 

it  is  necessary  and  sufficient  that  - (a^  sin  ¥0)  he  greater  than  the 

numerical  value  of  the  real  part  of  the  square  root  of  the  discriminant. 

For  this  to  be  true,  Marcos  - aj^sin  ¥Q)  must  be  negative. 

2 

shows  that  if  ifrQ  * “ 7 " x * (a^cos  ¥0  ~ /0)  “ (*3  + 

while  if  ¥ Q * - ♦ X -n  the  negative  of  this  appears.  Hence  stability  is 

possible  only  in  the  latter  case,  and  the  necessary  and  sufficient  conditions 
are: 

(21)  -f_  m - 7j  ♦ X « a,  a1  + a2  sin  ^ Q < 0. 

If  7/  » -n/2  approximately,  -f" 0 m ^ “ «/2 j or  if  X is  small,  the  nose 
is  yawed  in  about  the  same  direction  as  the  tail  is  bent. 


Substitution 

2 2. 1,2 
\ -*5> 
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An  oversimplified  physical  picture  is  given  in  Figure  where 
I2  ■ (y  ■ v'  - 0,  and  the  normal  force  on  the  tail  is  assumed  to  dominate 

the  others*  Then  7)  • - n/2  approximately.  The  projectile  is  vie.*)d  from 
the  rear  along  the  trajectory;  A,  T and  P are  the  projections  of  the  nose, 
the  trajectory  and  the  center  of  pressure  of  the  tail  on  the  plane  of  the 
paper.  Directed  arcs  represent  moments. 

Figure  5 showsthat,  in  the  stable  case,  an  increase  in  % causes  an 
opposing  increase  in  the  moment  due  to  the  normal  farce,  while  in  the 
unstable  case  this  moment  increases  Z still  more. 


A.  S.  GAIBRAITH 
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LIST  OF  CYM30LS 


Underlined  letters  represent  vectors;  e.g.,  u is  a vector  of  length  u. 

Aj  axial  moment  of  inertia. 

a^i  coefficients  in  "locking-in,,  equation;  defined  after  equation  (18). 

B:  transverse  moment  of  inertia. 

C,  D,  E,  F:  combinations  of  aerodynamic  symbols,  etc.  See  equations  (15) 

and  (17). 

d:  maximum  diameter  of  projectile. 

2 

Q:  (g/v  d)  sin  6 ♦ cos  6);  cross  component  of  gravitational  acceleration, 

essentially. 

2 

Gi*  °i  “ gid/u  ' 

g:  vector  acceleration  of  gravity;  £ ■ g-jX^  ♦ ^2^2  + 63*3  ♦ 
hj  distance  from  axis  to  center  of  pressure  of  tail,  in  calibers. 

Jj  j - pd^K/m. 

K:  dimensionless  aerodynanic  coefficient. 

KA~spin  damping  moment, 

K --moment  due  to  non-coincidence  of  lines  of  zero  lift  and  moment  of  tail, 
KDA--axial  drag, 

Kp— 3pin  accelerating  moment  due  to  canted  fins, 

Kp— Magnus  force, 

Ky— damping  (pitching)  moment, 

KL— lift, 

K^— static  moment, 

— normal  force, 

K^— cross-spin  (pitching)  force, 

K^— Magnus  moment, 

— Magnus  force  due  to  cross-spin. 
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LIST  CF  SYJBOLS  (Oon‘t) 

Mtgnus  moment  dv-c  to  cross-spin. 

Subscripts  1 and  2 refer  to  body  and  tail  respectively.  An  extra 
subscript  M is  a temporary  notation  to  indicate  the  moment  of  a force. 

2 2 

ki  transverse  radius  of  gyration,  in  calibers;  k ■ B/md  . 
k^j  axial  radius  of  gyration,  in  calibers. 

L>  distance  of  representative  point  from  c.g. 

l^i  direction  cosines  of  line  of  zero  normal  force  of  tail. 

iqi  mass  of  projectile. 

N:  (nutational)  yawing  rate;  N ■ (Am^/B)  (l+c). 

n.  s moment  caused  by  non-coincidence  of  lines  of  zero  lift  and  moment  of 
tail  acts  along  n^  ♦ ♦ n^* 

p:  cross  component  of  angular  velocity;  p ■ n exp  (i  co^da). 

Q:  rig)it  membsr  of  yawing  equation  (15),  when  gravity  is  neglected,  i3 


Q exp(i 


•J'  co^da); 


3ee  equation  (17). 


o..  ,q_:  terms  in  equations  of  motion  (lit*)  and  (llib)  respectively  that 

are  due  to  asymmetry. 

ra 

rs  a solution  of  the  yawing  equation  (17)  is  sought  xn  the  form  exp(  I rd  a). 

2 2 2 2 o 

3:  ballistic  (gyroscopic)  stability  f~clcr;  s = A k /1.E  J m. 

ti  distance  along  axis  from  c.g.  to  center  of  pressure  of  tail,  in  calibers. 

uj  vector  velocity  of  projectile;  u » u^x^  + UgXg  ♦ u^x^. 

V!  v * u/d,  Vj  ■ u^/d. 

axes  parallel  to  fixed  directions  in  space,  with  origin  at  c.g. 
x±t  unit  vectors  along  axes  x^,  which  are  fixed  in  projectile, 
z:  unit  vector  along  trajectory;  u ■ uz. 

z' j unit  vector  along  line  of  zero  lift  cf^tail, 
a:  distance  travelled  in  calibers;  a.  ■ J vdt. 


LIST  ls1  SYMBOLS  (Con’t) 


P:  angle  ol  misalignment  of  tail;  p » n/t. 


F • Og  ^ i . 

Ys  g/^2d. 

5j  angle  of  yaw,  from  trajectory  to  axis. 
£ 2 angle  of  c»nt  of  fins. 


£ : complex  yaw;  £ * exp(i 


ca^ da) , 


u g 2 x/^  2 2 

Tj  : angle  defined  by  sin  Tj  - a^/(a^  + a^  ) , cos  ij  ■ a^/ ( a^  + a^  ) , 

6:  angle  from  horizontal  to  trajectory,  positive  up. 

p 

Xs  approximately  measures  total  damping;  X ■ JN(l-A/b)  + (JH  + Jt)/k  . 
p,?  cross  angular  velocity;  p * 

2 2 2 2 , 

% 2 angle  defined  by  sin  X - a^/(a^  ♦ ) , cos  X - (*»^'  ♦ - a^  ) / 

(*3  * %*>  * 

^ 2 complex  yaw;  4 - - (v?  + iv^)/v. 

p:  density  of  the  air. 

2 , 

02  a ■ 1 - !/»• 


V2 


T ; snail  perturbation  of  "y"  • 

0:  Sulerian  angle  of  orientation. 
t(r  : Eulerian  angle  of  spin. 

f}  j angular  velocity  vector;  Qm  + ^2^2  + & j£.y 

“i*  Cl±/v* 
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